Spin-dependent electron-hole capture kinetics in conjugated polymers 
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The recombination of electron-hole pairs injected in extended conjugated systems is modeled as a 
multi-pathway vibron-driven relaxation in monoexcited state-space. The computed triplet-to-singlet 
ratio of exciton formation times r = tt/ts increases from 0.9 for a model dimer to 2.5 for a 32-unit 
chain, in excellent agreement with experiments. Therewith we rationalize recombination efficiency 
in terms of spin-dependent interstate vibronic coupling and spin- and conjugation-length-dependent 
exciton binding energies. 

PACS numbers: 78.60.Fi 73.61. Ph 82.20Wt 
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Introduction - The efficiency of electroluminescence 
(EL)from conjugated polymers in light emitting diodes 
(LEDs) is largely determined by the fraction of injected 
electron-hole (e-h) pairs which combine to form emissive 
spin-singlet (S) as opposed to non-emissive spin-triplet 
(T) bound states. @ If one assumes that the cross sec- 
tions for S and T e-h captures are equivalent as = or, 
then the singlet generation fraction \s = °~sl {o~s + 3<7t) 
will be limited to 25% according to spin multiplicity. 
Nonetheless, efficiencies corresponding to xs over 50% 
have been achieved in poly(p-phenylenevinylene) (PPV) 
based LEDs.||, || From this it has been inferred that 
singlet e-h capture is intrinsically more efficient than the 
respective triplet process, as > ot- Recent PA/PADMR 
experiments on a wide variety of conjugated systems indi- 
cate that the cross-section ratio r = as /<jt typically lies 
between 2 and 5, M and that the ratio increases with 
conjugation length. |5| 

The observed variation of r was originally related^) to 
the optical gap according to an interchain recombination 
model|^|; however, more recent experiments indicate that 
large values of r are characteristic of extended intrachain 
conjugation, j^, [7] Wilson et al. have measured xs and Xt 
consistent with r ~ 4 for a polymer, but only r w 0.9 for 
the monomer. While time-dependent scattering calcu- 
lations by Kobrak and Bittnerf^ indicate that the weakly 
bound S± exciton is a more efficient trap of free e-h pairs 
than the tightly bound Ti triplet exciton, the systematic 
variation of r with (effective) conjugation length n has 
remained unclear. 

This Letter addresses the kinetics of e-h capture in ex- 
tended conjugated systems. To our best knowledge, ours 
is the first molecular-based approach which consistently 
reproduces the spin- and conjugation-length dependence 
of the process, and rationalizes it in terms of exciton 
binding energies. This has important ramifications in 
the design and synthesis of polymer materials for device 
applications. 

Method - We simulate the recombination process in 
nondegenerate 7r-chains by the dissipative dynamics of 
the multi-level electronic system coupled to the phonon 
bath. The detailed structure of the Hamiltionian with a 
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is described in a separate publication,^] and we men- 
tion only those salient features that are relevant herein. 
The electronic part H e i reflects the configuration in- 
teraction (CI) of single excitations for a generic two- 
band polymer in a localized basis. (To[ [TTJ The configu- 
rations |m) = | mm) taken over valence and conduction- 
band Wannier functions |m) and \m) represent geminate 
(to = to) and charge-transfer (to ^ to) e-h pairs. hatH e i 
is parameterized from the 7r-band structure of extended 
PPV. Since the conjugated backbone is alternant, one- 
body CI integrals F mn = <W/m-n - $ m nf Wl -n of tnc 
band structure operator / obey electron-hole symmetry 
with f r = —f r . Spin dependency originates from two- 
body interactions with V mll = — (mn\\nm) + 2(mn\\mn) 
for singlets and V mll = —(mn\\nm) for triplets. As- 
suming zero-differential overlap except for geminate or- 
bitals, we take into account only true Coulomb and ex- 
change terms (toto||toto) and (mrn\\mm), and (transi- 
tion) dipole-dipole interactions (mn\ \mn) between gem- 
inate singlet pairs. The electron-phonon coupling term 
Hei-ph assumes that localized conduction/ valence lev- 
els ±/ Q and nearest-neighbor transfer integrals ±/i are 
modulated by lattice oscillations. Correspondingly, the 
phonon term consists of two sets of interacting local 
oscillators, giving rise to two dispersed optical phonon 
branches centered at 1600 and 100cm -1 , respectively. 
These frequencies roughly correspond to the C=C bond 
stretches and ring torsions in PPV, which largely dom- 
inate the Franck-Condon activity of the lowest optical 
transitions. |^2|, [l3|] Thus, by empirical adjustment of el- 
ph coupling strength in H e i- p h, we compute in Condon 
approximation |r| fairly accurate absorption and emis- 
sion band-shapes for PPV chains, reflecting conjugation- 
length dependence, vibronic structure and Stokes shifts. 

Diagonalizing H e i and H p h yields vertical excited 
states with energies e°|14 and normal modes with fre- 
quencies uj£, which allow us to transform the el-ph cou- 
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FIG. 1: Density of spin-singlet (left) and spin-triplet (right) 
states for PPV32 (shaded area) and their populations at differ- 
ent recombination stages. The t = population corresponds 
to the initially injected e-h pair. 

pling term in its diabatic representation. 

h = £^>>< a l + ^E(^i + p ?) 

+ E^Qda)(b|. (2) 

FIG.l shows the S and T densities of states (DOS) for a 
conjugated chain with n = 32 repeat units (PPV32), ob- 
tained as a convolution of Lorentzian line-shapes with 
r = O.OleV. We note that in singlet state-space, op- 
tical coupling with the ground state is concentrated 
mainly around the DOS bottom since only S± and other 
low-lying excitonic states contain geminate configura- 
tions with transition dipoles in additive combinations. |o| 
Nonequilibrium dynamics in excited state-space is de- 
rived further from the diabatic interstate vibronic cou- 
plings g° b £- Assuming that phonons thermalize rapidly 
on the time scale of electronic relaxation, internal conver- 
sion (IC) rates k a \, for one-phonon processes are obtained 
in Markov approximation. 

(g° ) 2 

kab = TtV" °^ fog + l)(r(Ct>£ - UJab) - r(cj£ + W a6 ))(3) 

Here nj is the Bose-Einstein distribution of phonons at 
T = 300K and T is the empirical broadening. Note 
that for an elementary a — ► b conversion to occur there 
must a phonon mode uj^ close to the transition frequency 
uj a i, = (e°—e%)/H. Singlet and triplet superposition states 
are separately evolved in time according to the quantum 
master equation for the density matrix p in energy rep- 
resentation, 

Pab = —i^abPab — E Rab,cdPcd (4) 
cd 

where we decouple populations p aa from coherences p a b 
according to the Bloch model. fig] 

RaaM = —k a b + S a b k ac (5) 

c 



Rab.ab = g ^(kac + k bc ) (6) 

c 

Results - Equations (3-6) allow us to model the relaxation 
of an arbirtary monoexcitation, including that with an 
electron and hole injected into the far ends of the conju- 
gated chain. The resulting initial population of the DOS 
is shown in FIG.l. We show as well the populations 
at an early relaxation time (40ps), and at 640ps when 
a stationary state is reached. In both spin-spaces, half 
of the population density propagates down to the lowest 
excitons (Sxt or Txt),[|16| but the other half remains 
locked in metastable charge-transfer (CT) states (Sct or 
Tqt)- The branching of the relaxation pathways is due 
to e-h symmetry, which separates excited states into even 
and odd under e-h transposition. The XT states, which 
are the lowest even ones, are not vibronically coupled 
with the CT states, which the lowest odd ones. XT and 
CT states are ultimately populated equally because the 
initial separated e-h pair is a 1:1 mixture of eigenstates 
with even and odd e-h parity. The weak e-h symmetry 
breaking of real conjugated systems is further modeled by 
slight conduction/ valence band nonsymmetry according 
to the relation f r /f r = —1.1. The perturbation barely 
affects the DOS, but invokes weak vibronic coupling be- 
tween the states of odd and even types, so that they 
undergo slow mutual conversions. The time-dependent 
XT and CT populations arc shown in FIG. 2 both for the 
symmetric (a) and nonsymmetric (b) cases. We see that 
when e-h symmetry is broken, the fast formation of XT 
and CT during the the first 100-200ps is followed up by 
a slow CT— »XT relaxation. Most notably, in comparison 
with the symmetric case, the initial XT-CT branching ra- 
tio changes drastically in favor of Sxt for singlets and in 
disfavor of Txt for triplets. The subsequent Sct —> Sxt 
relaxation is also substantially faster than the respective 
Tct — > Txt conversion process, and Sxt population re- 
mains about twice as high as that of Txt throughout 
the simulated time-range. We note here that CT states 
are barely spin-dependent (see FIG.l) due to negligible 
spin-exchange; hence S ^=> T intersystem crossing pre- 
equilibrium of long-lived CT states is also likely to occur 
prior to final e-h binding. As shown in FIG. 2a, Sxt for- 
mation proceeds notably faster than that of Txt- Since 
even under e-h symmetry population buildup deviates 
from first-order kinetics in early stages when intermedi- 
ate states are formed (see FIG.l), we measure the rel- 
ative efficiency of exciton formation by the time r at 
which XT population reaches 40%. Thus for PPV32 we 
find ts — 109ps and tt = 263ps, which corresponds to 
r = tt/ts = 2.5 or xs — 45% in good agreement with 
both EL|, § and PA/PADMR|, § data. Singlet ver- 
sus triplet recombination enhancement is related to the 
difference in exciton binding energies, |^| which we mea- 
sure by e b = e° CT — s XT . Whereas Sxt is slightly lower 
than Sct, spin-exchange contributes ~leV more to the 
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FIG. 2: Essential-state populations as a function of time for 
PPV32 with strict (a) and weakly broken (b) e-h symmetry. 
The Sct and Tct populations not shown in (a) follow closely 
the Sxt curve. 
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FIG. 3: (a) Variation of singlet and triplet binding ener- 
gies and formation times with inverse conjugation length 1/n. 
(ts — tt ~ o, e b s = ■, e g = □). (b) Binding energy and 
formation-time ratio tt/ts • and binding-energy ratio £s/st 
■ vs. 1/n. Inset: r vs. e s /ey. Dashed curves are guides for 
the eye. 



triplet binding energy e\. Therefore relaxation into Txt 
takes on the average a longer sequence of internal con- 
versions, whose energy cutoff at «0.2eV is determined by 
the phonon spectrum. If we assume that relaxation time 
is proportional to the dissipated energy and reciprocal to 
the mutual coupling of the states along the pathway, we 
can write down, 



CX- ' vib h 



(7) 



where r V ib reflects any difference of the effective vibronic 
coupling in S and T state-space. We further probe re- 
lation (7) by varying the conjugation length n. For the 
sake of brevity, results for the series of PPV-model chains 



with n = 2, 4, 8, 16, and 32 are summarized in FIG. 3a, 
where exciton binding energies e b and formation times 
r are plotted vs. inverse chain length 1/n as commonly 
adopted. Note that e b s and e T decrease in parallel with 
1/n and converge yet for the longest chains, suggest- 
ing that effective conjugation length is over 10 repeat 
units. Exciton formation times show different behavior. 
Whereas tt increases about four times on going from the 
dimer to PPV32, ts increases only slightly. In FIG. 3b 
we plot r vs. 1/n and establish an excellent agreement 
with experiments^], [?| for both oligomers and polymers, 
including a notable triplet enhancement with r — 0.9 for 
n = 2. The inset is a plot of r vs. e T /e b s , which illustrates 
the existence of a rough linear relationship (7) with r v n, 
similar to 0.6. 

Conclusion - In summary, we have proposed and mod- 
eled an intrachain c-h recombination mechanism for con- 
jugated systems, with computational results consistent 
with the spin-resolved experiments in all essential points. 
The model reflects the notion that exciton formation 
occurs via vibron-mediated internal conversions across 
excited states, and hence the overall rate of electronic 
relaxation is inversely proportional to the energy dissi- 
pated into vibrons. Hence while strong vibronic cou- 
pling of the lowest triplet is responsible for r < 1 in 
short chains and molecular species, the disproportion of 
exciton binding energies accounts for the large r values 
in extended polymers. The approximate e-h symmetry 
of conjugated systems implies the parallel formation of 
long-lived charge-transfer states, which are nearly spin- 
independent, and thus susceptible to intersystem crossing 
prior to final spin-singlet enhanced e-h binding. 
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